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We present a systematic study of charged pion and kaon interferometry in Au+Au collisions 
at ^Sjy„=200 GeV. The kaon mean source radii are found to be larger than pion radii in the 
outward and longitudinal directions for the same transverse mass; this difference increases for more 
central collisions. The azimuthal-angle dependence of the radii was measured with respect to the 
second-order event plane and similar oscillations of the source radii were found for pions and kaons. 
Hydrodynamic models qualitatively describe the similar oscillations of the mean source radii for 
pions and kaons, but they do not fully describe the transverse-mass dependence of the oscillations. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

Measurements of the quark-gluon plasma (QGP), pro¬ 
duced in nucleus-nucleus collisions at the Relativistic 
Heavy Ion Collider ^HIC) [H-I^ and the Large Hadron 
Collider (LHC) [l-Q, showed that the QCP exhibits 
rapid hydrodynamic expansion, followed by hadroniza- 
tion, which results in the emission of many particles. The 
time of last scattering among hadrons is referred to as ki¬ 
netic freeze out. To understand the dynamics and prop- 
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erties of the QCP, it is important to understand the full 
system evolution and how it is constrained by the mea¬ 
surements of the space-time distribution at kinetic freeze 
out. 

The quantum statistical interferometry of identical 
particles, also known as Hanbury Brown Twiss (HBT) 
interferometry or femtoscopy, is a powerful tool to mea¬ 
sure the spatial and temporal scales of systems cre¬ 
ated in nucleus-nucleus collisions i i. This tech¬ 
nique was first developed to measure the angular di¬ 
ameter of stars through intensity interferometry of ra¬ 
dio waves 0- It has also been applied to nuclear and 
particle physics m- In nucleus-nucleus collisions, the 
interferometry using emitted hadrons measures the spa¬ 
tial extent of the particle-emitting source at the time of 
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kinetic freeze out. 

Despite the successful description of various observ¬ 
ables at RHIC by the hydrodynamic models P, Q , there 
was significant discrepancy between HBT data and the¬ 
oretical models P, [l^ ■ Recent theoretical development 
has improved the agreement by including realistic physics 
conditions such as stiffer equation of state and a viscosity 
of the created matter [l^ . 

Charged pions are often used for the interferometry 
analysis because of their abundant production, but re¬ 
cently acquired large data sets by RHIC and LHC experi¬ 
ments allow study of the particle-species dependence [H- 
Kaon interferometry is of particular interest, because 
the contribution from resonance decays is reduced com¬ 
pared to pions mill, thereby providing a more direct 
view of the particle-emitting source. PHENIX at RHIC 
published an analysis of one-dimensional source imaging 
for charged kaons M- STAR at RHIC has recently 
published three-dimensional source imaging [l^ . where 
charged kaons lack the nonCaussian tail in the source 
function observed in the pion sample. This result may 
be caused by the reduced contribution from long-lived 
resonances as well as a different time dependence due to 
a shorter rescattering phase. Further systematic stud¬ 
ies using different particle species are needed to better 
constrain the space-time evolution and freeze-out distri¬ 
butions of the created medium. 

The HBT measurement is also sensitive to the initial 
spatial anisotropy and the subsequent evolution of the 
created matter. Due to the strong collective expansion, 
one might expect the eccentricity of the source shape in 
the initial state to be reduced at freeze out and possibly 
to be reversed if the collective expansion is stronger in the 
direction of the reaction plane, or if the expansion time is 
sufficiently long. To probe the spatial source anisotropy 
at freeze out, HBT measurements with respect to the 
event planes have been performed using two-pion corre¬ 
lation [iqI - I^ . Large oscillations of the pion source radii 
relative to the second-order event plane were observed, 
which indicates that the pion source at freeze out is elon¬ 
gated in the direction perpendicular to the event plane 
even after the collective expansion. 

In this paper, we present azimuthal-integrated and 
azimuthal-dependent source radii for charged pions and 
kaons in Au-|-Au collisions at ^/s^= 200 GeV. Results 
are compared with the hydrodynamic models for both 
particle species. 


II. EXPERIMENT 

The PHENIX experiment [2^ is designed to mea¬ 
sure particles produced in nucleus-nucleus collisions with 
good momentum resolution, including photons, electrons, 
muons, and hadrons, to study properties of the QGP. The 
PHENIX detectors are comprised of magnet systems and 
detectors for particle tracking and identification, event 
timing, plus vertex position and centrality determina¬ 


tion. The particle tracking and identification detectors 
are arranged into central and forward (muon) arms. Fig¬ 
ure [T] shows the layout of the PHENIX detector during 



FIG. 1. (Color online) The layout of PHENIX detectors in the 
2007 run configuration. The top figure shows the central arm 
detectors viewed along the beam axis. The bottom figure 
shows the side view of the global detectors and muon arm 
detectors. 

Global detectors characterize the global event charac¬ 
teristics in heavy ion collisions. The beam-beam counters 
(BBG) measure the collision time, and the position 
of the collision vertex along the beam axis, as well as the 
collision centrality. The BBC comprises two identical sets 
of counters located ±144 cm from the nominal collision 
point and surrounds the beam pipe covering the pseu¬ 
dorapidity range of 3.0 < \r]\ < 3.9. Each BBC has 64 
modules of Cerenkov radiators, and measures the number 
of charged particles in its acceptance. The zero-degree 
calorimeters are located 18 m from the nominal collision 
point and measure the energy of s pec tator neutrons. The 
reaction-plane detector (RXNP) was installed prior 
to the 2007 RHIC run to measure the event-plane angle in 
heavy ion collisions. The RXNP comprises two sets of 24 
scintillators on both north and south sides, and is located 
±39 cm from the vertex position. The scintillators are 
arranged around the beam pipe in two concentric rings 
of 12 segments in the azimuthal direction. The outer and 
inner rings cover pseudorapidity ranges of 1.0 < \r]\ < 1.5 
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and 1.5 < |? 7 | < 2.8, respectively. 

The PHENIX central arms comprise two sets of detec¬ 
tors located on the west and east sides of the beam axis. 
Each arm covers 90° in azimuth and a pseudorapidity 
range of I 77 I < 0.35. Track and momentum reconstruc¬ 
tions of charged particles were performed with the drift 
chambers (DC) and pad chambers (PC). The DC are lo¬ 
cated at a radial distance of 2.02 m to 2.46 m from the 
beam axis in the west and east arms, covering 2 m length 
along the beam axis. The PC are multi-wire proportional 
chambers in each of the central arms, and are located at 
radial distances of 2.5 m (PCI) and 4.9 m (PCS). The 
tracks and momenta were reconstructed by combining 
the hit information in the DC and PCI, providing a mo¬ 
mentum resolution of 6p/p ~ 1.3%0l.2% xp GeV/c [ 2 ^ . 
Global-track reconstruction was performed by associat¬ 
ing these tracks with hits in the outer detectors, such 
as PCS and the lead-scintillator (PbSc) electromagnetic 
calorimeters, as shown in Fig.[T] Particle identification is 
provided by the PbSc which is a sampling calorime- 
ter with a timing resolution of about 500 ps located 
at a radial distance of 5.1 m from the beam axis. 


III. DATA ANALYSIS 


The squared mass of the particles is given by the follow¬ 
ing formula: 



where p is the momentum, t is the time-of-flight, L is 
the flight path length, and c is the speed of light. Pi- 
ons and kaons were selected from a 2cr window around 
their peaks in the squared mass distribution. Additional 
requirements, i.e. to be away from the mass peak of 
other particles, were applied to reduce contamination. 
The tt/A separation was achieved up to a momentum 
of ^1 GeV/c. Contamination in the pion samples from 
kaons is below 1% for p r; 1 GeV/c and contamination in 
the kaon samples from pions (protons) is below 4% (1%) 
for p Ri 1 GeV/c. 

C. Construction of the Correlation Function 

In this section, a bold character denotes four¬ 
dimensional vector and an arrow denotes three- 
dimensional vector. 

The experimental correlation function defined as 


The ^/s^ = 200 GeV Au-I-Au collision data were col¬ 
lected by PHENIX during the 2007 running period. A 
total of 4.2 billion events were used for this analysis, 
where the minimum bias trigger with at least two hits in 
each BBG was required. This trigger measures 92±3% 
of the total inelastic cross section . Additional offline 
requirements of one zero-degree-calorimeter hit on each 
side and a collision vertex position of less than ±30 cm 
were applied. 


A. Track Selection 

Charged tracks with good quality were selected based 
on the track information from the DC and PCI. To re¬ 
duce the background due to the random association of 
hits and reconstructed tracks, track residuals were re¬ 
quired to be less than 2a in the cj) — z plane at the PC3 
and PbSc for pions. For kaons, this cut was relaxed to 
2.5(7 to increase statistics. The fraction of the random 
background is ~4.6% after the 2a cuts and ~5.3% af¬ 
ter the 2.5(7 cuts at pt =0.5 GeV/c in the 0%-10% most 
central collisions. The effect of the track quality cuts was 
included in the systematic uncertainty. 


B. Particle Identification 

Particle identification was performed by combining 
time-of-flight data from the PbSc in the west arm, the 
reconstructed momentum, and flight path length from 
the collision vertex to the hit position at the PbSc wall. 


C^2(q) 


^(q) 

BiciV 


( 2 ) 


was measured as a function of the pair momentum differ¬ 
ence q = Pi — P2 , where A and B are constructed from 
identical particle pairs from the same event and mixed- 
event respectively. The mixed-events are taken from sim¬ 
ilar event centralities and vertex positions. In the case of 
azimuthal-dependent analysis, the mixed-events are also 
required to have similar values for the second-order event 
plane defined in Sec lIIIDl 

Particle pairs with similar momenta and spatially close 
to each other are affected by incorrect track reconstruc¬ 
tion and detector inefficiencies. These effects were re¬ 
moved by applying pair selection cuts at the DC and 
PbSc following our previous analysis [13) 113 ■ In addition, 
pairs that are associated with hits on the same tower of 
the PbSc were removed. 

The particle pairs were analyzed with the Bertsch- 
Pratt parameterization [30l| as functions of the pair 
momentum difference q and mean pair momentum fc, 
where k = {pi ±p^)/2. The k is projected into its longitu¬ 
dinal component k^ and transverse component fcx. The q 
is projected into the longitudinal (qi), outward (go), and 
sidewards {qs) components, where qi denotes the beam 
direction, go is perpendicular to kx, and qs is perpendic¬ 
ular to both gi and go. In this frame, the energy (tem¬ 
poral) component of the four-dimensional vector is taken 
in the outward component by performing the analysis in 
the longitudinal co-moving system, where kz = 0. 

The (72 (q) function is divided into two components 
based on the core-halo picture in which the A parameter 
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controls the relative strength of the core and the halo. 
C72(q) = 

= X{l + G{q))F,iq) + {l-X) (3) 

G{qs, go, gi) = (4) 

The Fc{q) is the Coulomb correction factor evaluated by 
the Coulomb wave function [^, , where g is the scalar 

quantity of q. The central core contributes to the quan¬ 
tum statistical interference. The halo includes the decay 
of long-lived particles for which the quantum statistical 
interference occurs in a g range that is too small to be re¬ 
solved experimentally, and for which the Coulomb inter¬ 
action is negligible. The core is assumed to be a Gaussian 
source as given by Eq. Q. 

The HBT radii denoted by Rs, i?o, and Ri represent 
the spatial extent of the emission region in each direc¬ 
tion, but Ro and to a lesser extent R\ and i?os include 
a contribution from the emission duration. All radii are 
sensitive to position-momentum correlations. The i?os 
term arises in the case of azimuthal-dependent analysis 
due to asymmetries in the emission region Q, while it 
vanishes in the azimuthal-integrated analysis. 

The HBT radii were measured as a function of /ct 
and presented as a function of the transverse mass 
rriT = (fcr^ + to study particle-species depen¬ 

dence, where m is the particle mass. 

D. Event Plane Dependence 

The second-order event-plane angle (4>2) was deter¬ 
mined using the RXNP detector based on the azimuthal 
anisotropy of emitted particles in momentum space: 

Qx = ^w,cos(2(t>i), (6) 

Qy = '^w^s\rL{24>i), (7) 

where (jji is the azimuthal angle of each segment i in the 
RXNP and Wi is the weight which reflects the particle 
multiplicity in that segment. Corrections for detector 
acceptance as detailed in Ref. [s^ were applied. 

Due to the finite number of particles within the RXNP 
acceptance, the observed event plane 4^2 is smeared 
around the true event plane $. This smearing effect 
is typically accounted for by the resolution. The event 
plane resolution defined as Res{4t2} = (co s[2( 4'2 — $)]) 
was estimated by the two-subevent method (3J| using the 
event plane correlation between the RXNP at forward 
and backward angles. The Res{4'2} has a maximum of 
0.75 in midcentral events [2^. 

The finite event plane resolution reduces the oscillation 
amplitude of HBT radii relative to the event plane. To 
take this effect into account, a mo del-independent correc¬ 
tion suggested in Ref. [s^ was applied. The pair distribu¬ 
tion measured at a certain azimuthal angle (j) relative to 


the reconstructed event plane, N{q,(j) — '£' 2 ), is smeared 
by the finite event plane resolution and finite width of 
angular bins A. The Fourier coefficients for the true and 
measured N{q,(j) — 4'2), Na,n{q, (t> ~ ^ 2 ), can be associ¬ 
ated with the following relation: 

jV-P(q, 0 - VI/ 2 ) = <-(q, (/.-$) ''''^2/2 (8) 

X (cos[n(4'2 - 4))]), 

where a denotes sine and cosine terms of the Fourier co¬ 
efficients {a = s,c) and n denotes the order of the coeffi¬ 
cient. The above equation is analogous to the correction 
for the elliptic flow (^ 2 ). Based on Eq. (jH]), the A(q) and 
B(q) functions can be unfolded by using the following 
equation: 

nbln/2 

-^(9,'5^j) “ -^exp(g,'5^j) + 2 y ^ 

n—7n,2m,--- 

X [^c,n^('7) cos{n(j)j) + iV®^(g) sm{n(j)j)], 

(9) 

where ribin is the number of azimuthal angular bins, and 
m is the order of the event plane, and (f>j denotes the 
center of angular bin which corresponds to the az¬ 
imuthal angle of the pair with respect to the event plane, 
and Cn.m(A) are given by 

= (^exp(g, (k - ^ 2 ) COs[n{(j) - ^- 2 )]), 

— y ^ A0xp(g, ) cos(yz(^j)/?rbin, 

j 

= (^exp(g, (/> - ^ 2 ) sin[n((/) - 4'2)]), 

= A^exp(g, (l>j) sm{n(j)j)/nhin, (H) 
j 

nAj2 

sin(nA/2)(cos[n('I'5^® — 4))]) 

The details of Eqs. (in))- (IT^ can be found in Ref. [s^. 


E. Systematic Uncertainties 

Systematic uncertainties were estimated by variations 
of track quality cuts at PCS and PbSc, pair selection 
cuts, and particle identification (PID) cuts. Also, the 
effect of the Coulomb correction was studied by varying 
the input source size in the calculation of Fc{q) in Eq. ([3|). 
Typical systematic uncertainties of the measured radii for 
charged pions and kaons are listed in TableUand Tablellll 
In the azimuthal-dependent analysis, the variations 
when using different event planes from forward, back¬ 
ward, and both combined RXNPs, were also incorpo¬ 
rated. The systematic uncertainties of the oscillation 
amplitudes of HBT radii were dominated by the event 
plane determination, which were 16% on average in the 
final eccentricity defined by the oscillation of Rs^ and the 
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TABLE I. Typical systematic uncertainties of HBT parame¬ 
ters for positive pion pairs in 0%-10% centrality and 0.6 < 
kx < 0.7 GeV/c. 

systematic source 

A 

R. 

Ro Ri 



[%] 

[%] [%] 

track quality 

1.8 

0.3 

0.5 3.1 

pair selection 

4.3 

1.0 

4.6 3.7 

Particle ID 

0.4 

0.3 

1.3 0.0 

Coulomb 

0.4 

0.1 

0.3 0.1 

Total 

4.7 

1.1 

4.8 4.8 


TABLE II. Typical systematic uncertainties of HBT parame¬ 
ters for charge combined kaon pairs in 0%-10% centrality and 
0.3 <kT < 0.68 GeV/c. 

systematic source 

A 

R. 

[%] 

Ro 

[%] 

Ri 

[%] 

track quality 

5.1 

2.2 

1.9 

2.2 

pair selection 

9.0 

1.5 

0.1 

1.9 

Particle ID 

6.1 

0.3 

4.5 

0.1 

Coulomb 

4.6 

0.3 

1.1 

0.2 

Total 

12.9 

2.7 

5.0 

2.9 


same fraction of the uncertainty was assumed for pions 
and kaons. 

The effect of momentum resolution was studied em¬ 
ploying the same method as previous analyses [H, [s^ ■ 
The momentum was smeared according to the known mo¬ 
mentum resolution and the correlation function was re¬ 
constructed using the smeared A(q) and -B(q) functions. 
By taking the ratio of the smeared and unsmeared corre¬ 
lation function, the correction factor was obtained. The 
correction on the momentum resolution was performed 
by multiplying the correction factor to the measured cor¬ 
relation function. The correction did not affect Rs and 
i?i, but slightly increased A (<10%) and Ro (<6%). 


IV. RESULTS AND DISCUSSION 
A. Azimuthal-integrated analysis 

Figure [2][a-c) shows an example of correlation func¬ 
tions of pion pairs and kaon pairs in 0%-10% centrality 
in a kr bin with fit lines given by Eq. where the 
momentum correction is not applied. The kx range is 
selected to have similar mx for pions and kaons. The 
3-dimensional >l(q) and B{q) functions are projected in 
each q direction. In the projection, the other q are re¬ 
stricted to be less than 40 MeV/c (e.g. when making 
C 2 ( 9 s), the projection ranges of qo and qi should be go 
<40 MeV and q\ <40 MeV). The ID correlation func¬ 
tions shown in Fig. [5] are obtained by taking ratio of the 
projected A(q) and i?(q) functions. The extracted HBT 


radii with the statistical uncertainties are also shown in 
each panel. The width of the enhancement at the low q 
region in the correlation function is proportional to the 
inverse of the HBT radius. The width of the correla¬ 
tion function is comparable between pions and kaons in 
the sidewards direction, but narrower for kaons in the 
outward and longitudinal directions, indicating a larger 
radius in those directions than for pions with a similar 
mean uit- We note that the data points at lower q bins 
fluctuate due to lower statistics. The effect of the fluc¬ 
tuation on the radii was studied by varying the fit range 
and it was found to be within a few percent for both 
pions and kaons. 


PHENIX Au+Au 200 GeV 

0-10% 

Rs=(3.49±0.04) fm 
^ Rs=(3.78±0.15)fm 

0 TiTi ,|3) 

(m))=0.79 GeV/c 
R;;=(3.02±0.04) fm 
Ro=(3.99±0.15) fm 

• K^KbKK’ (j.) 

(ml)>=0.76 GeV/c 
R|‘=(2.95±0.04) fm 
rI^=( 4.07±0.18) fm 

B 





0 0.1 0 0.1 0 0.1 

q [GeV/c] q [GeV/c] q [GeV/c] 

S 0 I 


FIG. 2. (Color online) Correlation functions of negative pi¬ 
ons and charged kaons for 0%-10% centrality (a-c), where 
positive and negative kaons are combined. Open boxes show 
the systematic uncertainties. Solid and dashed lines show the 
fit functions and the extracted radii values are shown in the 
figure. 

Figure |3] shows the extracted HBT parameters of 
charged pions and kaons for four centrality classes as a 
function of rriT- Results for charged pions in the low 
rriT region from STAR [s^ are also plotted. The source 
parameters from the two experiments are in good agree¬ 
ment, but the A parameters are 20% lower at low mx. 
The value of A is sensitive to the combinatorial back¬ 
ground level, which may differ between PHENIX and 
STAR. Positive and negative pions are quite consistent. 
The presented data are also consistent with our previous 
results MM- 

The decrease of HBT radii with mx is often attributed 
to the position-momentum correlation induced by col¬ 
lective flow. The slope of the mx dependence becomes 
steeper for more central collisions, which is consistent 
with an expectation of a stronger radial flow Q . Rs shows 
approximate mx scaling between pio ns and kaons as pre¬ 
dicted by the Buda-Lund model [4^, which is based on 
the analytic approach of the perfect fluid hydrodynamics. 
On the other hand, Ro and R\ of kaons show larger values 
than those of pions as noted already in Pig. [U where the 
mx scaling is broken. The difference increases with cen¬ 
trality going from peripheral to central collisions. The 
similar difference between pions and kaons for i?i was 
reported by STAR [T^. 

The results are compared with the hydrokinetic model 



















FIG. 3. (Color online) Extracted HBT parameters of charged pions and kaons as a function of ttit for the centralities indicated, 
where open boxes show the systematic uncertainties. Results of charged pions from STAR are compared. Calculations 
from the hydrokinetic model (HKM) and viscous-hydrodynamic model (Bozek) are also shown. 


(HKM) [3i,|4ll. The HKM incorporates realistic con¬ 
ditions such as the Glauber initial condition, crossover 
transition, fluid hydrodynamics, microscopic transport 
and resonance decays, but does not explicitly include the 
viscous correction. It is reported that the model calcu¬ 
lations with the initial condition of the color glass con¬ 
densate are very similar to those with the Glauber initial 
condition [dlj. As shown in Fig. [3l the HKM de¬ 
scribes well the overall trend of HBT radii for pions and 
kaons in all centrality bins, however it overestimates Ro 
of pions in more central collisions and underestimates Rs 
and R\ of pions in peripheral collisions. The HKM also 
describes the difference of pions and kaons in the lon¬ 
gitudinal direction, which can be understood by strong 
transverse flow , but the difference in the outward di¬ 
rection cannot be explained well. The data for pions in 


most central collisions are also compared with (3-|-l)-D 
viscous hydrodynamic model [s^ calculations which em¬ 
ploy a Glauber initial condition and 77 /s = 0.08 (also see 
Sec. IIV B 31 for details). The model follows the general 
trends in the data. 

The ratio of Ro and i?s, which is sensitive to the emis¬ 
sion duration of particles, is also plotted as a function 
of rriT in Fig. |4l Results for both species do not show 
any significant centrality dependence, but the values for 
kaons are larger than those for pions at all ttit and cen¬ 
tralities, a possible indication of longer emission duration 
time for kaons than for pions. The HKM reproduces the 
data for kaons well, but not for pions. 

The ttit scaling of HBT radii was inspired by the hy¬ 
drodynamic expansion (^ . This is based on the idea 
that the kinetic freeze out of hadrons occurs at the same 
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FIG. 4. (Color online) Ratio of and Rg for charged pions and kaons as a function of niT- Calculations from the hydrokinetic 
model (HKM) [s^ and viscous-hydrodynamic model (Bozek) are also shown. 



FIG. 5. (Color online) HBT radii of charged pions and kaons as a function of kr, where open boxes show the systematic 
uncertainties. 


time and the hadrons with similar velocities are emitted 
from the same homogeneity region. In other words, the 
homogeneity length depends on the particle mass under 
the presence of radial flow. In Fig.jSJ both pion and kaon 
HBT radii for central and peripheral events are plotted as 
a function of kx- Unlike the case of the rriT dependence 
shown in Figl^l both radii seem to be scaled better for kr 
in all q directions as predicted in Ref. [d^. This model 
includes many different effects such as the hadronic cas¬ 
cade and resonance decays in addition to radial flow. 

We have also checked charge-dependent kaon HBT 
radii in Fig. [51 There was no significant difference be¬ 
tween positive and negative kaons as we expected. If nu¬ 
cleons are dominant in the particle-emitting source and 
the net baryon density is not small, the measured radii 
might be different between K'^ and K~ (and also pi¬ 
ons) because of smaller cross section of K~^ — N than 


K — N [i^. However this is not observed. 


B. Azimuthal-dependent analysis 

1. Results 

We have measured the azimuthal angle dependence of 
HBT radii with respect to 'I '2 for both charged pions 
and kaons. Figure [Zl[a-c) shows the correlation func¬ 
tions of charged kaons in the 20%-60% centrality bin 
in the in-plane (|(/) — 'I'2| < 7r/16) and out-of-plane 
(1^ — 'I'2 — 7r/2| < 7r/16) directions without correction 
for the event plane resolution. The correlation functions 
in Fig. [7] are calculated in the same way as Fig. [U i.e. 
when making the 1-dimensional C 2 along the q of inter¬ 
est, the other q are limited to be less than 50 MeV/c. 
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FIG. 6. (Color online) Comparison of HBT radii between 
positive and negative kaons pairs in central and peripheral 
collisions, where open boxes show the systematic uncertain¬ 
ties. 



FIG. 7. (Golor online) Correlation functions of charged kaons 
in 20%-60% centrality (a-c), where positive and negative 
kaons are combined. The correlation functions along qs and 
Qo directions are averaged out between positive and negative 
q. Lines show the fit functions and the extracted radii values 
are shown in the figure. 


To make a comparison of the C 2 width between in-plane 
and out-of-plane directions, the C 2 in the positive and 
negative q are averaged because they are symmetric over 
q = 0 within the statistical uncertainties. The extracted 
radii without the correction are also shown in the figure. 
A difference of the width in the correlation function be¬ 
tween these (in- and out-of-plane) directions can be seen 
in the sidewards (Fig.[7Ka)) direction. Figure [5] shows the 
extracted HBT radii of charged kaons as a function of az¬ 
imuthal pair angle (j) with respect to 'F 2 for two centrality 
bins where (kr) is ^0.77 GeV/c. We first fix A in Eq. ([3]) 
by taking the average for A obtained in all azimuthal bins, 
then we fit in individual azimuthal bins again with fixed 
A parameter as detailed in [s^ . This treatment is based 


FIG. 8. (Color online) Squared HBT radii of charged kaon 
pairs as a function of azimuthal pair angle (j> with respect to 
5'2 for two centrality bins, where fcr is integrated over 0.3-1.5 
GeV/c. The open symbols at <() —’F 2 = tt are the same data as 
that at (^ — ^2 = 0. Open boxes show systematic uncertainties 
and the solid lines are the fitting functions given by Eq. m- 


on the assumption that A has no azimuthal angle depen¬ 
dence, and the data fluctuate but do not depend on the 
azimuthal angle beyond the systematic uncertainty. The 
cosine oscillations of and Ro^ (Fig. |SKa,b)) and the 
sine oscillation of i?os^ (Fig- Eld)) can be clearly seen. 
Non-zero i?os at {(j) — '£' 2 ) = f implies that the di¬ 
rection of the particle emission is tilted relative to the 
main axis of the emission region. The oscillation of Rg^ 
seems to be larger than for i?o^ in 20%-60% centrality 
bin. 

We have also measured the charged pion HBT radii 
with respect to 'I '2 for the same centrality bins as kaons 
with six kx {ttit) bins as shown in Fig. The averages 
of Rg^, Ro^, and decrease with kx as seen in Fig. [31 
The Rg^ and Ro^ have similar but opposite oscillations 
in all kr bins. For 20%-60% centrality, both transverse 
radii show finite oscillation even in the lowest kx bin, 
which indicates that the pion emission happens from an 
elliptical source. For 0%-20% centrality, the Rg^ has a 
weak azimuthal angle dependence, while the Ro^ has a 
larger oscillation than the Rg^. It could be consistent 
with i?o^ being more influenced by the anisotropic flow 
as discussed in our previous publication The oscil¬ 
lations of Rog^ decrease with kx and Ri displays a neg¬ 
ligible azimuthal angle dependence, which qualitatively 
agree with hydrodynamic calculations [ssl . [dg . 

The data shown in Fig.jSjand Fig. [3] are fitted with the 
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FIG. 9. (Color online) Squared HBT radii of charged pion pairs as a function of azimuthal pair angle (j) with respect to ^2 for 
6 fcr bins and 2 centrality bins ((a)-(d) for 0%-20% and (e)-(h) for 20%-60%), where open symbols at 0 — ^'2 = tt is the same 
data as that at <)) — ^'2 = 0. Open boxes show systematic uncertainties and solid lines show the fit functions by Eq. (II3. 


functions below to extract the oscillation strength (d^ : 

+ 2 i ?^,2 cos{2/\(t>) = s, o, 1), 

i?2(A</))=2i?2 2sin(2A</)) (fi = os), (13) 

where 2 the 2''‘^-order Fourier coefficient and 
A(j) = (/) — 4'2- Detailed discussion on the oscillation 
amplitudes is presented in Sec. IIV B 31 


2. Blast-wave model fit 

In this section, we perform blast-wave model fits to our 
results to extract features at the kinetic freeze out and 
study their particle species dependence. The blast-wave 
(BW) model [d^ is based on a hydrodynamical model 
parameterized by the freeze-out conditions, such as the 
freeze-out temperature (T/) and the transverse flow ra¬ 
pidity (po)- This model is further expanded in Ref. [s^ 
to describe the elliptic flow and azimuthal angle depen¬ 
dence of HBT radii by introducing additional parameters: 
2nd-order modulation in transverse flow rapidity (p 2 ), the 
transverse source size (i?x, Ry), the freeze-out time (tq), 
and the emission duration (Ar). Once the above seven 
parameters are fixed, the pt spectra, elliptic flow, and 
HBT radii can be calculated within the model. 

Each freeze-out parameter has a different sensitivity 
to each experimental observable (35j |. For example, the 
P 2 and the ratio of i?x and Ry are less sensitive to pt 


spectra, but more sensitive to the elliptic flow and the 
azimuthal angle dependence of HBT radii. To effectively 
constrain those parameters, a fit to pt spectra was first 
performed to determine T/ and po, then the other param¬ 
eters were determined by simultaneous fit to the elliptic 
flow and azimuthal-dependent HBT radii. For the source 
size parameters, the i?x and Ry/Rx are actually used as 
the fitting parameters. 

The BW model assumes that the freeze out for all 
hadron species takes place at the same time, but the ac¬ 
tual situation may be more complicated. To investigate 
how the extracted freeze-out parameters vary by particle 
species of HBT radii, the following fits were tested: 

A. Fit for Pt spectra and V2 of tt, K, p along with 
HBT radii of tt. 

B. Fit for Pt spectra and V2 of tt, K, p along with 
HBT radii of K. 

In the case of Fit B, both azimuthal-dependent and 
azimuthal-integrated HBT radii of charged kaons were 
included in the fit. 

Figure [TPl shows the results of Fit A for the pt spectra 
(Fig. fTUT a)) and elliptic flow (Fig. fTUI d)) of tt, K and p, 
and pion HBT radii (Fig. fTUf b.c.e.f)) in 20%-60% cen¬ 
trality. The solid lines show the BW fitting functions and 
the dashed lines in the panel (a) and (d) represent the 
extended fitting function beyond their actual fit ranges. 
Only three kT bins for HBT radii are shown here, but all 
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FIG. 10. (Color online) The blast-wave model fits (Fit A) to the pr spectra (a), elliptic flow of tt, K, p (d), and HBT radii of 
TV (b,c,e,f). The solid lines show the fit functions. 


six kx bins shown in Fig. Owere simultaneously used in 
the fit. Here, the data of pr spectra are taken from 
and the data of V 2 from [49|. The px spectra and V 2 are 
well described at low px, and the overall trend of the 
HBT radii is also reproduced by the BW model. 

The results from the fits are summarized in Table imi 
The systematic uncertainties of the BW model fit were 
estimated by varying the fit conditions: the fit range, 
a surface diffuseness to control the density profile 
and the relative weighting factor between different par¬ 
ticle species. The systematic uncertainties of data were 
taken into account in the calculation of Our results 
from Fit A are in good agreement with those in previous 
studies [ 33 , 113 . The results of Fit B shows slightly dif¬ 
ferent values, i.e, smaller and larger At. The smaller 
source (i?x) might be intuitively understood as due to 
kaons freezing out earlier than pions, but it is not signif¬ 
icant in the parameter r. The At obtained by the Fit B 
using the kaon HBT result shows relatively larger values 
than the results by Fit A, which is consistent with the 
result from Ro/Rs as shown in Fig. |4l 

Also, the Tot dependence of the pion and kaon HBT 
radii has been calculated using the parameters obtained 
from Fit A as shown with lines in Fig. 1111 For a com¬ 
parison, the 0*^-order Fourier coefficients (i?/i,o) for pi- 
ons which correspond to the HBT radii obtained in the 


azimuthal-integrated analysis are plotted as filled sym¬ 
bols. The kaon HBT radii from the azimuthal-integrated 
analysis are also compared in the figure. The BW model 
shows the tv jK difference in the sidewards and outward 
directions, but not in the longitudinal direction unlike 
the experimental data. 


3. Oscillation amplitudes with hydrodynamic models 

The BW model [s^ suggests that the source eccen¬ 
tricity at freeze out is given by Efinai = 2 i?g 2 /^s 0 = 
—2i?Q 2 /-Rs,o = 2 i?Qg_ 2 /^s.o foP- dSl) fo the absence 
of position-momentum correlation, i.e. radial flow. In 
the presence of radial flow, the above relation would be 
smeared because the HBT radius does not reflect the 
whole source size, but the esnai from R^ 2 could still be a 
good estimator in the limit oi kx = 0. 

The HBT radii of both pions and kaon averaged over 
the azimuthal direction are well described by the hydro- 
dynamic models including the BW model as shown so far. 
In this section, the oscillation amplitudes are also com¬ 
pared with the hydrodynamic models. The oscillation 
amplitudes were extracted by using Eq. m- The sys¬ 
tematic uncertainties were estimated by performing the 
fitting with Eq. dn for the data of various systematic 
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TABLE III. Summary of extracted parameters in the blast-wave model fit for two fitting conditions (see the text for details). 
The Tf and po parameters were obtained by fits to pr spectra, and the other parameters were obtained by a simultaneous fit 
to V 2 and the HBT radii. The values in parentheses represent the systematic uncertainties derived by varying the model fit 
conditions. 


Fit 

Centrality 

Tf 

[MeV] 

po 

P2 

Ax 

[fm] 

Ay/Ax 

[fm/c] 

Ar 

[fm/c] 

xVndf 

(spectra) 

xVndf 

{V2) 

xVndf 

(HBT) 

A 

0%-20% 

104 

0.995 

0.047 

11.28 

1.092 

8.22 

2.06 

143.4/27=5.3 

25.2/21=1.2 

2526.7/96=26.3 



(5) 

(0.055) 

(0.005) 

(0.23) 

(0.003) 

(0.23) 

(0.18) 




A 

20%-60% 

113 

0.905 

0.074 

8.25 

1.171 

6.1 

1.56 

206.5/27=7.6 

46.4/21=2.2 

1998.3/96=20.8 



(8) 

(0.059) 

(0.012) 

(0.29) 

(0.003) 

(0.31) 

(0.12) 




B 

0%-20% 

104 

0.995 

0.042 

10.19 

1.102 

8.48 

2.73 

143.4/27=5.3 

11.0/21=0.5 

116.9/28=4.2 



(5) 

(0.055) 

(0.004) 

(0.46) 

(0.004) 

(0.81) 

(0.58) 




B 

20%-60% 

113 

0.905 

0.067 

7.44 

1.182 

4.95 

2.86 

206.5/27=7.6 

39.4/21=1.9 

95.7/28=3.4 



(8) 

(0.059) 

(0.01) 

(0.55) 

(0.008) 

(0.91) 

(0.52) 





E 



E 


o 

or 



171^- [GeV/c] 



m^- [GeV/c] 
PHENIX Au+Au 200 GeV 
71 0-20% Q 71 20-60% 

K 0-20% K 20-60% 

BW 71 0-20%- - BW 7t 20-60% 

BW K 0-20% BW K 20-60% 


m-r [GeV/c] 



FIG. 11. (Color online) The 0‘''-order Fourier coefficient of 
HBT radii of charged pions as a function of tut, and blast- 
wave model calculations for both pions (dashed line) and 
kaons (solid line), where the blast-wave model parameters 
shown in Fit A of Table [Till were used. 


pendence especially in 0%-20% centrality. It should be 
noted that the Ro and i?os contain the particle emission 
duration in addition to the geometrical information, and 
they are also dominated by the anisotropy in the expan¬ 
sion velocity. 

The data are compared with the blast-wave calcula¬ 
tion using the Fit A parameters in Fig. [121 The de¬ 
pendency on rriT of the oscillation amplitudes is not de¬ 
scribed well, although the rriT dependence of the mean 
radii is reproduced well. The large of HBT in Fit 
A in Table [ml is mainly caused by such a discrepancy. 
The calculations from the event-by-event 3-l-l-D viscous- 
hydrodynamic model [ 3 ^ with a Glauber initial condition 
and shear viscosity ? 7 /s= 0.08 (and also nonzero bulk vis¬ 
cosity) are compared to the same data in Fig. [131 The 
model employs the equation of state with a crossover 
transition, but does not include a hadron cascade. The 
model quantitatively agrees with the data of Rs and i?os 
but overestimates Ro in 20%-60% although the trend 
(not the magnitude) of dependency on tot is reproduced 
in contrast with the blast-wave model. 

When data points of kaons are compared to those of pi¬ 
ons in the 0%-20% centrality sample, they mostly agree 
within the systematic uncertainties, but in the 20%-60% 
sample the 2i?s 2 /-Rs,o( 2 -Ro. 2 /-^o,o) kaons is slightly 
larger(smaller) than that of pions. More precise mea¬ 
surement is needed to confirm the difference. 


sources described in Sec lIII El In Fig. [12] and FigjTH the 
oscillation amplitudes with 4 different combinations of 
HBT radii are plotted in the form of a final eccentricity, 
2i?^ 2 /-^y,o> where p and v denote o,s,os. The i?^ 2 is 
a fitting parameter in Eq. (USD that can take a negative 
value, which represents a different phase of the cosine 
function and shown in Fig. fT^ b.c) and Fig. fT^ b.c). 

Figure USDa) shows the oscillation amplitude of i?s^ 
relative to the average, which is most sensitive to the fi¬ 
nal source eccentricity. The value of 2R‘^ 2 /^s,o increases 
with nriT, which would reflect tot dependent elliptic- 
ity of the emission region. The other combinations of 
|2i?^ 2 /-R^, 2 l show similar tot dependence, but less de- 


V. SUMMARY AND CONCLUSION 

We have presented results from the PHENIX exper¬ 
iment on charged pion and kaon femtoscopy measure¬ 
ments in Au-I-Au collisions at = 200 GeV. In 

the azimuthal-integrated analysis, we have measured the 
HBT radii of both species with fine tot and centrality 
bins. Tot scaling holds well for R^, but there are visi¬ 
ble differences for Ro and Ri between charged pions and 
kaons at the same tot, and the differences become larger 
in more central collisions, tot scaling breaks for those 
radii, but kp scaling works well for all radii. It is ob- 
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FIG. 12. (Color online) Oscillation amplitudes relative to the average for four different combinations of the azimuthal-dependent 
HBT radii as a function of rtiT for charged pions and kaons. Open boxes show systematic uncertainties. Calculations from the 
blast-wave model with parameters of Fit A shown in Table ITIIl are shown for comparison. 



FIG. 13. (Color online) Oscillation amplitudes relative to the average for four different combinations of the azimuthal-dependent 
HBT radii as a function of mr for charged pions. Open boxes show systematic uncertainties. Calculations from the 3-1-TD 
viscous-hydrodynamic model are shown for comparison. 


served that the ratio Ro/Rs of kaons is larger than that 
of pions, which may imply different emission durations. 
The hydrokinetic model was compared with our data. 
It reproduces most aspects of the data of both charged 
pions and kaons, but it fails to accurately describe the 
difference in /Jo¬ 
in the azimuthal-dependent analysis, a first measure¬ 
ment of the HBT radii of charged kaons with respect 
to the second-order event plane has been performed and 
compared with pion measurements with finer rriT bins. 
Oscillation with respect to the event plane of kaon HBT 
radii has been clearly observed, and is similar to that 
of pions. The data were compared with the blast-wave 
model and the 3-l-l-D viscous-hydrodynamic model. The 
blast-wave model provides a good description of the over¬ 
all trend of the pr spectra, the elliptic flow, and the 
mean HBT radii, but fails to describe the details of fem- 


toscopy measurements, such as the mx dependent os¬ 
cillation amplitude of the source radii. While the 3-1-1- 
D viscous-hydrodynamic model does qualitatively repro¬ 
duce the data, it overestimates the oscillation of Ro- We 
note that the viscous hydrodynamic model also repro¬ 
duces well the other observables such as the px spectra 
and elliptic flow [5l|. 

Both the hydrokinetic model and viscous- 
hydrodynamic model surprisingly describe all aspects 
of the femtoscopic data, even though these models 
lack the shear viscosity of plasma and the microscopic 
transport phase. Including these effects may improve 
the description of the measured Ro/Rs and /?o, 2 - More 
precise measurements and systematic model comparison 
for both azimuthal-dependent and azimuthal-integrated 
HBT measurements are needed. The particle-species 
dependence, in addition to the differential femtoscopy 
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measurements may help to elucidate the expansion 
dynamics of heavy ion collisions. 
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Appendix A: Blast-wave model 

An emission function in the blast-wave parameteriza¬ 
tion for bosons [s^ is given by 

S{r, (j)s,T, T]) = tut cosh(77 - y)D(r, 

OO 

X ^nacos{<j>t-4>p) cosh{r]-y) 

n—1 

(Al) 


defined as 

Oi = ^smhp{r,(j)s), (A2) 

fi = ^coshp{r,(j)s). (A3) 

The transverse rapidity p(r, (fig) is defined as: 

p(r, (/)s) = (po + P 2 cos(2(^f,))f, 

f = ^(cr/i?.)2 + {y/Ry)^, (A4) 

where x and y are space coordinate of particles, /ps is 
azimuthal angle of the spatial positions, and (j)b is a boost 
direction. It is assumed that particles are boosted to 
the direction perpendicular to the elliptical subshell of 
the particle-emitting source, which satisfies the relation 
below: 


tan((/is) = (Ry/Rx)^tan((/)b). (A5) 

The distribution of the source elements D(f) is given by 
D(f) = l/(l + e("-i)/“), (A6) 


where a denotes a surface diffuseness and a = 0 gives a 
box profile and a = 0.3 gives approximately a Gaussian 
profile. Observables, such as spectra, V 2 , and HBT radii, 
are obtained by performing the integral of the emission 
function Eq. m over phase space weighted with certain 
quantity B: 


j d'^x, S{x, K)B{x, K) = 


(A7) 





rdrS^r, 


(f>g,T, T])B{x,K). 


Azimuthally integrated pr spectra can be obtained by 
integrating over (j)p and r in Eq. (IA8I) setting B{x, K)=l. 
If we assume Boltzmann distribution for all particles, 
only the first term in the summation in Eq. mil is used. 
Also, in case of analyzing particles in midrapidity region, 
Eq. (I All) can be simplified by setting y=0. Then Eq. (IA8I) 
can be rewritten as the following: 


dN 

PxdpT 


= y/^ToAr 


d(l)p 


d(j)s 


rdr 


dp 


X TOTC0sh(7?)D(r,(/),)e““"('^‘-'^^)e-^™"^(’’), 


= TqAt 


r27T 


r27T 


d(J}p 


d(j)s 


rdr 




(A8) 


where Kn{[3) is the modified Bessel function of the second 
kind, which is defined as 


Kn{z) 


1 

-J dtcosh(nt)e-^™"‘'(‘^ 


(A9) 


where uit is the transverse mass, y is rapidity, (j)p is az- Here we replace (j)b — (j)p as (j)’, and the range of the inte- 
imuthal angle of particle momentum, and a and P are gral over (j)’ is from pb to pb — 27r. Then the range can 
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be replaced from 0 to 27r because the integrand is the 
periodic function with Finally, Eq. (IA8I) is rewritten 
as 


dN 

prdpT 


2(27r)^/^roArmT 


X 




rdr n{r, 4>s) Io{a)Ki{l3), 


(AlO) 


where In is the modified Bessel function of the first kind 
given by 


In{z) = 


27r 


dt cos(nt)e 


— z cos(t) 


(All) 


The elliptic flow V 2 is calculated as 


f d4>p f d'^x cos{2(j)p)S{x,K) 


(A12) 


The denominator is the same expression with Eq. (lAlOl) . 
The numerator can be calculated by a similar way to 
derive the pr spectra. 


d(j)p J d^x cos{2(j)p)S{x, K), 


p27r n27r poo 

= 2V^toAt / d(j)p / d(j)s / rdr 0s) 

do do do 

p27r poo p2TT 

= 2\/^roAr / d0p / rdr mT ^^(r, 0s)^i(/5) cos( 205 ) / d0'cos(20/)e"^°®^'^ \ 

do do do 

p27T poo 

= 2{2tt)^^‘^ toAttut / d(j)p / rdr r2(r, 0s)Ari(/3) cos( 2 (/)b)/ 2 (a). 

Jo Jo 

Finally, the elliptic flow can be expressed as 

/ \ jQ''d(j)pJ^rdrft{r,(j)s)Ki{P)cos{2(j)b)l2{a) 

V2(pT,rn) = -—-. 

lo rdrn{r,(l)s)Io{a)Ki{l3) 

The HBT radii are related to space-time variance as [s^ 


+ (y^)) “ “ {f))cos{2(j)p) - (iy)sin(20p), 

+ if)) + “ (d^)) cos{2(j)p) + (iy) sin(20p), 

-2(3T{{ix) cos(j)p + {iy) sin^p) -|- 


1 


Ros = {xy) cos{2(t)p) - ) - (y ))sin(2(/)p) +/3T((ti) sinc^ip - (ty) cos((>p). 


2 
s 

= {z^)-2Pi{tz)+0f{P), 

= (5"), 


2 /T2\ 


(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 


where 


(/(^)) = 


/ d‘^xf{x)S{x, K) 
f d‘^xS(x, K) 


= x>^-{x^), 


(A19) 

(A20) 


Appendix B: Data tables 

The extracted HBT radii and the oscillation ampli¬ 
tudes for charged pion and kaons in Au-|-Au collisions at 
y/s^ =200 GeV are summarized in Table IIVIIXII 


and Pi vanishes in the LCMS frame and the terms includ¬ 
ing t and z depend on the proper time r and emission 
duration of particles At. As shown in above equations, 
Rs depends on only the spatial extent of the source and 
azimuthal angle pp, while Rq and i?os are also sensitive 
to T and At as well as the spatial extent. 













17 


TABLE IV. HBT parameters of positive pion pairs, shown as value ± statistical uncertainty [absolute value] ± systematic 
uncertainty [%] for the centrality bins shown in Fig. [S] 


Centrality 

(mr) 

[GeV/c] 

A 

R, 

[fm] 

Ro 

[fm] 

Ri 

[fm] 

0%-10% 

0.3 

0.292 ± 0.003 ± 11.1 

5.11 ± 0.03 ± 0.5 

5.55 ± 0.04 ± 2.9 

6.3 ± 0.04 ± 2.3 


0.36 

0.358 ± 0.004 ± 8.7 

4.76 ± 0.03 ± 0.6 

5.22 ± 0.03 ± 3 

5.41 ± 0.04 ± 2.8 


0.41 

0.392 ± 0.004 ± 6.1 

4.52 ± 0.03 ± 0.8 

4.91 ± 0.03 ± 2.8 

4.8 ± 0.04 ± 3.7 


0.47 

0.402 ± 0.005 ± 11 

4.28 ± 0.03 ± 1.3 

4.51 ± 0.03 ± 4.7 

4.31 ± 0.03 ± 5.2 


0.53 

0.45 ± 0.006 ± 7.3 

4.17 ± 0.03 ± 0.7 

4.21 ± 0.04 ± 4.2 

4.08 ± 0.04 ± 4.2 


0.59 

0.454 ± 0.007 ± 5.9 

3.98 ± 0.03 ± 1.4 

3.97 ± 0.04 ± 5.7 

3.8 ± 0.04 ± 3.2 


0.66 

0.458 ± 0.007 ± 4.7 

3.77 ± 0.03 ± 1.1 

3.62 ± 0.04 ± 4.8 

3.45 ± 0.04 ± 4.8 


0.79 

0.462 ± 0.008 ± 6.1 

3.47 ± 0.03 ± 1.2 

3.27 ± 0.04 ± 4.9 

3.01 ± 0.04 ± 5.2 


0.96 

0.447 ± 0.016 ± 12 

3.16 ± 0.06 ± 3.1 

2.91 ± 0.07 ± 10.9 

2.49 ± 0.06 ± 6.1 


1.16 

0.448 ± 0.034 ± 18.5 

2.93 ± 0.11 ± 8.3 

2.4 ± 0.12 ± 14.4 

2.13 ± 0.11 ± 9.4 

10%-20% 

0.3 

0.313 ± 0.003 ± 9 

4.63 ± 0.03 ± 0.6 

5.01 ± 0.03 ± 2.4 

5.77 ± 0.04 ± 1.9 


0.36 

0.379 ± 0.004 ± 7.2 

4.3 ± 0.03 ± 0.5 

4.74 ± 0.03 ± 2.6 

4.95 ± 0.04 ± 2.8 


0.41 

0.41 ± 0.004 ± 7.3 

4.1 ± 0.03 ± 1.3 

4.41 ± 0.03 ± 3.1 

4.44 ± 0.03 ± 3.5 


0.47 

0.441 ± 0.005 ± 7.5 

3.94 ± 0.03 ± 1 

4.12 ± 0.03 ± 3.6 

4.07 ± 0.03 ± 4.2 


0.53 

0.487 ± 0.006 ± 6.6 

3.83 ± 0.03 ± 0.8 

3.88 ± 0.03 ± 2.9 

3.83 ± 0.03 ± 3.7 


0.59 

0.501 ± 0.008 ± 5.5 

3.67 ± 0.03 ± 0.9 

3.68 ± 0.04 ± 3.7 

3.56 ± 0.04 ± 2.3 


0.66 

0.501 ± 0.008 ± 4.2 

3.5 ± 0.03 ± 1.5 

3.38 ± 0.03 ± 3.7 

3.24 ± 0.03 ± 3.5 


0.79 

0.501 ± 0.008 ± 4.3 

3.21 ± 0.03 ± 0.9 

2.99 ± 0.03 ± 4.9 

2.78 ± 0.03 ± 3.1 


0.96 

0.515 ± 0.016 ± 8.2 

2.94 ± 0.05 ± 1.4 

2.78 ± 0.06 ± 4.1 

2.34 ± 0.05 ± 6.3 


1.16 

0.52 ± 0.035 ± 11.2 

2.76 ± 0.1 ± 9 

2.23 ± 0.1 ± 2.8 

2.09 ± 0.1 ± 9.7 

20%-40% 

0.3 

0.339 ± 0.003 ± 7.2 

4.09 ± 0.02 ± 0.6 

4.34 ± 0.02 ± 2.5 

5.01 ± 0.03 ± 1.5 


0.36 

0.401 ± 0.003 ± 5.8 

3.81 ± 0.02 ± 0.9 

4.09 ± 0.02 ± 2 

4.35 ± 0.03 ± 2.1 


0.41 

0.433 ± 0.004 ± 5.4 

3.63 ± 0.02 ± 0.7 

3.87 ± 0.02 ± 2.3 

3.88 ± 0.03 ± 2.2 


0.47 

0.46 ± 0.005 ± 6.9 

3.5 ± 0.02 ± 0.9 

3.59 ± 0.02 ± 3.2 

3.55 ± 0.03 ± 3.7 


0.53 

0.493 ± 0.005 ± 5.8 

3.39 ± 0.02 ± 0.8 

3.38 ± 0.02 ± 2.9 

3.3 ± 0.03 ± 3 


0.59 

0.519 ± 0.007 ± 3.7 

3.26 ± 0.02 ± 0.3 

3.2 ± 0.03 ± 2.7 

3.11 ± 0.03 ± 2.1 


0.66 

0.505 ± 0.006 ± 3.2 

3.07 ± 0.02 ± 0.4 

2.91 ± 0.02 ± 2.7 

2.77 ± 0.03 ± 1.9 


0.79 

0.532 ± 0.007 ± 3.5 

2.89 ± 0.02 ± 1.1 

2.66 ± 0.02 ± 3.7 

2.45 ± 0.02 ± 3.1 


0.96 

0.54 ± 0.014 ± 7.7 

2.63 ± 0.04 ± 4.3 

2.34 ± 0.04 ± 3.5 

2.07 ± 0.04 ± 3.2 


1.16 

0.554 ± 0.027 ± 4.3 

2.35 ± 0.08 ± 2.9 

2.09 ± 0.08 ± 5.5 

1.76 ± 0.07 ± 5.6 

40%-70% 

0.3 

0.363 ± 0.004 ± 5.8 

3.27 ± 0.02 ± 0.9 

3.39 ± 0.03 ± 2.6 

3.94 ± 0.03 ± 1.1 


0.36 

0.426 ± 0.005 ± 2.8 

3.05 ± 0.03 ± 0.6 

3.29 ± 0.03 ± 1.7 

3.42 ± 0.03 ± 1.1 


0.41 

0.455 ± 0.006 ± 4.2 

2.92 ± 0.03 ± 0.7 

3.16 ± 0.03 ± 3 

3.09 ± 0.03 ± 1.5 


0.47 

0.493 ± 0.007 ± 3.6 

2.87 ± 0.03 ± 1.2 

2.94 ± 0.03 ± 1.7 

2.89 ± 0.03 ± 2.1 


0.53 

0.509 ± 0.008 ± 3.9 

2.73 ± 0.03 ± 0.9 

2.81 ± 0.03 ± 1.8 

2.63 ± 0.03 ± 2 


0.59 

0.54 ± 0.01 ± 3.4 

2.7 ± 0.03 ± 1.2 

2.67 ± 0.04 ± 1.7 

2.46 ± 0.04 ± 1.6 


0.66 

0.569 ± 0.01 ± 1.2 

2.61 ± 0.03 ± 2.1 

2.48 ± 0.03 ± 1.9 

2.31 ± 0.03 ± 1 


0.79 

0.551 ± 0.01 ± 4.6 

2.36 ± 0.03 ± 1.6 

2.1 ± 0.03 ± 2.1 

1.93 ± 0.03 ± 2.1 


0.96 

0.619 ± 0.024 ± 5.7 

2.35 ± 0.06 ± 1 

1.99 ± 0.06 ± 4.4 

1.82 ± 0.06 ± 1.6 


1.16 

0.592 ± 0.041 ± 3.8 

2.02 ± 0.13 ± 2.8 

1.64 ± 0.13 ± 11.1 

1.45 ± 0.11 ± 3.5 
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TABLE V. HBT parameters of negative pion pairs, shown as value ± statistical uncertainty [absolute value] ± systematic 
uncertainty [%] for the centrality bins shown in Fig. [S] 


Centrality 

{mr) 

[GeV/c] 

A 

Rs 

[fm] 

Ro 

[fm] 

Ri 

[fm] 

0%-10% 

0.3 

0.275 ± 0.003 ± 10.5 

5.11 ± 0.03 ± 0.6 

5.49 ± 0.04 ± 3 

6.28 ± 0.05 ± 2.9 


0.36 

0.353 ± 0.004 ± 7.3 

4.74 ± 0.03 ± 0.7 

5.25 ± 0.04 ± 2.1 

5.44 ± 0.05 ± 3.4 


0.41 

0.387 ± 0.005 ± 7.9 

4.5 ± 0.03 ± 0.4 

4.91 ± 0.04 ± 3.6 

4.82 ± 0.04 ± 3.9 


0.47 

0.399 ± 0.005 ± 9.1 

4.31 ± 0.03 ± 0.9 

4.38 ± 0.04 ± 4.3 

4.29 ± 0.04 ± 5.4 


0.53 

0.444 ± 0.007 ± 4.7 

4.17 ± 0.03 ± 1.2 

4.19 ± 0.04 ± 3.6 

4.11 ± 0.04 ± 3 


0.59 

0.45 ± 0.008 ± 7.4 

3.95 ± 0.04 ± 1.2 

3.91 ± 0.04 ± 5.7 

3.82 ± 0.04 ± 3.8 


0.66 

0.451 ± 0.008 ± 7.3 

3.77 ± 0.03 ± 0.2 

3.61 ± 0.04 ± 5.3 

3.47 ± 0.04 ± 3.3 


0.79 

0.442 ± 0.008 ± 6.2 

3.49 ± 0.03 ± 1 

3.16 ± 0.04 ± 6 

2.96 ± 0.04 ± 5.5 


0.96 

0.437 ± 0.017 ± 5 

3.22 ± 0.07 ± 2.4 

2.78 ± 0.07 ± 4.5 

2.46 ± 0.06 ± 6.9 


1.15 

0.526 ± 0.046 ± 13.3 

3.1 ± 0.13 ± 3.1 

2.27 ± 0.13 ± 8.1 

2.43 ± 0.13 ± 11.1 

10%-20% 

0.3 

0.295 ± 0.003 ± 9.4 

4.63 ± 0.03 ± 0.4 

4.93 ± 0.04 ± 3.4 

5.74 ± 0.04 ± 2.3 


0.36 

0.363 ± 0.004 ± 6.4 

4.3 ± 0.03 ± 0.3 

4.71 ± 0.04 ± 3.1 

4.88 ± 0.04 ± 2.6 


0.41 

0.404 ± 0.005 ± 6.9 

4.13 ± 0.03 ± 0.5 

4.42 ± 0.04 ± 2.8 

4.46 ± 0.04 ± 3.6 


0.47 

0.443 ± 0.006 ± 7.2 

3.95 ± 0.03 ± 1 

4.13 ± 0.03 ± 2.9 

4.05 ± 0.04 ± 4.1 


0.53 

0.462 ± 0.007 ± 5.9 

3.8 ± 0.03 ± 0.2 

3.91 ± 0.04 ± 3.4 

3.77 ± 0.04 ± 2.8 


0.59 

0.492 ± 0.008 ± 4.1 

3.68 ± 0.03 ± 0.4 

3.64 ± 0.04 ± 3.7 

3.53 ± 0.04 ± 3.5 


0.66 

0.48 ± 0.008 ± 4.7 

3.47 ± 0.03 ± 1.2 

3.33 ± 0.03 ± 4.3 

3.15 ± 0.03 ± 3.8 


0.79 

0.497 ± 0.009 ± 5.1 

3.25 ± 0.03 ± 0.5 

3.01 ± 0.03 ± 4 

2.75 ± 0.03 ± 3.7 


0.96 

0.482 ± 0.016 ± 3.6 

2.89 ± 0.05 ± 1.2 

2.68 ± 0.06 ± 4.5 

2.27 ± 0.05 ± 5.2 


1.15 

0.579 ± 0.043 ± 16.5 

2.67 ± 0.12 ± 8.6 

2.29 ± 0.11 ± 5.7 

2.11 ± 0.1 ± 6.7 

20%-40% 

0.3 

0.32 ± 0.003 ± 6.8 

4.11 ± 0.02 ± 0.5 

4.28 ± 0.03 ± 2.3 

5.05 ± 0.03 ± 2 


0.36 

0.385 ± 0.004 ± 7.3 

3.83 ± 0.03 ± 1.4 

4.03 ± 0.03 ± 2.8 

4.32 ± 0.03 ± 2.1 


0.41 

0.429 ± 0.005 ± 5 

3.65 ± 0.02 ± 0.5 

3.89 ± 0.03 ± 1.3 

3.93 ± 0.03 ± 2.5 


0.47 

0.458 ± 0.005 ± 6.6 

3.5 ± 0.02 ± 0.5 

3.59 ± 0.03 ± 3.1 

3.54 ± 0.03 ± 3.5 


0.53 

0.487 ± 0.006 ± 4.5 

3.4 ± 0.03 ± 0.7 

3.38 ± 0.03 ± 1.8 

3.32 ± 0.03 ± 2.7 


0.59 

0.509 ± 0.007 ± 3.5 

3.26 ± 0.03 ± 0.9 

3.18 ± 0.03 ± 2.7 

3.1 ± 0.03 ± 2.9 


0.66 

0.521 ± 0.007 ± 2.5 

3.11 ± 0.02 ± 0.6 

2.96 ± 0.03 ± 2.6 

2.86 ± 0.03 ± 2.2 


0.79 

0.521 ± 0.007 ± 4.8 

2.87 ± 0.02 ± 0.7 

2.62 ± 0.03 ± 3 

2.4 ± 0.02 ± 3 


0.96 

0.536 ± 0.015 ± 2.2 

2.66 ± 0.04 ± 1.1 

2.31 ± 0.04 ± 3 

2.05 ± 0.04 ± 3.3 


1.16 

0.565 ± 0.034 ± 8.6 

2.45 ± 0.1 ± 3.4 

1.97 ± 0.09 ± 1.7 

1.81 ± 0.09 ± 3.3 

40%-70% 

0.3 

0.344 ± 0.004 ± 5.2 

3.25 ± 0.03 ± 0.9 

3.31 ± 0.03 ± 1.5 

3.96 ± 0.04 ± 0.7 


0.36 

0.409 ± 0.006 ± 4.3 

3.05 ± 0.03 ± 0.9 

3.27 ± 0.03 ± 1.9 

3.41 ± 0.04 ± 1.8 


0.41 

0.452 ± 0.007 ± 3.1 

2.92 ± 0.03 ± 1 

3.15 ± 0.03 ± 2.3 

3.1 ± 0.04 ± 1.3 


0.47 

0.489 ± 0.008 ± 4.8 

2.83 ± 0.03 ± 1.4 

2.95 ± 0.03 ± 2.6 

2.86 ± 0.04 ± 1.8 


0.53 

0.528 ± 0.009 ± 2.5 

2.8 ± 0.03 ± 1.1 

2.77 ± 0.03 ± 2 

2.67 ± 0.04 ± 2.2 


0.59 

0.539 ± 0.011 ± 2.9 

2.65 ± 0.04 ± 2.3 

2.62 ± 0.04 ± 3.2 

2.47 ± 0.04 ± 1.3 


0.66 

0.532 ± 0.01 ± 4.6 

2.5 ± 0.03 ± 2 

2.4 ± 0.03 ± 2.6 

2.25 ± 0.03 ± 3.7 


0.79 

0.565 ± 0.012 ± 3.5 

2.41 ± 0.03 ± 1.4 

2.19 ± 0.04 ± 1.7 

1.99 ± 0.03 ± 2.5 


0.96 

0.582 ± 0.023 ± 7.7 

2.26 ± 0.07 ± 2.2 

1.91 ± 0.07 ± 3.9 

1.7 ± 0.06 ± 3.5 


1.15 

0.641 ± 0.063 ± 3.5 

2.18 ± 0.17 ± 9 

1.72 ± 0.14 ± 10.9 

1.67 ± 0.15 ± 17.3 
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TABLE VI. HBT parameters of charge-combined kaon pairs, shown as value ± statistical uncertainty [absolute value] ± 
systematic uncertainty [%] for the centrality bins shown in Fig. |3l 


Centrality 

{mr) 

[GeV/c] 

A 

R. 

[fm] 

Ro 

[fm] 

Ri 

[fm] 

0%-10% 

0.76 

0.507 ± 0.044 ± 12.9 

3.8 ± 0.15 ± 2.7 

4.07 ± 0.15 ± 5 

4.09 ± 0.17 ± 2.9 


0.91 

0.46 ± 0.043 ± 12.7 

3.32 ± 0.13 ± 4.9 

3.49 ± 0.15 ± 5.3 

3.25 ± 0.17 ± 4.1 


1.1 

0.404 ± 0.043 ± 4.3 

2.86 ± 0.14 ± 3.4 

3.1 ± 0.18 ± 5.4 

2.36 ± 0.15 ± 7.7 

10%-20% 

0.76 

0.539 ± 0.044 ± 11.3 

3.37 ± 0.13 ± 2.6 

3.56 ± 0.13 ± 4 

3.55 ± 0.16 ± 4.5 


0.91 

0.515 ± 0.044 ± 15.3 

3.03 ± 0.12 ± 2.7 

3.2 ± 0.13 ± 1.2 

2.98 ± 0.14 ± 12.5 


1.1 

0.457 ± 0.052 ± 14.7 

2.75 ± 0.14 ± 1.5 

2.48 ± 0.14 ± 2.6 

2.34 ± 0.17 ± 10.8 

20%-40% 

0.76 

0.64 ± 0.044 ± 11.4 

3.18 ± 0.11 ± 2 

3.28 ± 0.1 ± 4 

3.15 ± 0.12 ± 4.7 


0.91 

0.489 ± 0.034 ± 7.6 

2.69 ± 0.1 ± 2.8 

2.73 ± 0.1 ± 3 

2.29 ± 0.11 ± 3.8 


1.09 

0.501 ± 0.046 ± 13.9 

2.56 ± 0.13 ± 3.8 

2.31 ± 0.12 ± 3.2 

2.05 ± 0.13 ± 7.4 

40%-70% 

0.76 

0.624 ± 0.062 ± 5.1 

2.73 ± 0.16 ± 5.3 

2.67 ± 0.16 ± 2.5 

2.35 ± 0.15 ± 8.3 


0.91 

0.565 ± 0.065 ± 14.6 

2.33 ± 0.17 ± 9.8 

2.44 ± 0.16 ± 5.5 

1.97 ± 0.17 ± 15.4 


1.09 

0.575 ± 0.068 ± 19.5 

2.08 ± 0.17 ± 5.8 

2.03 ± 0.19 ± 5.3 

1.54 ± 0.16 ± 8 


TABLE VII. Ro/Rs of positive and negative pion pairs plus charge-combined kaon pairs, shown as value ± statistical uncertainty 
[absolute value] ± systematic uncertainty [%] for the centrality bins shown in Fig. [d] 


meson 

pair 

(mr) 

[GeV/c] 

0%-10% 

Ro j R‘. 

10%-20% 

20%-40% 

40%-70% 

7r"'"7r+ 

0.3 

1.09 ± 0.01 ± 3.3 

1.08 ± 0.01 ± 2.3 

1.06 ± 0.01 ± 2.1 

1.04 ± 0.01 ± 1.9 


0.36 

1.10 ± 0.01 ± 3.5 

1.10 ± 0.01 ± 3.0 

1.07 ± 0.01 ± 1.6 

1.08 ± 0.01 ± 1.7 


0.41 

1.09 ± 0.01 ± 3.5 

1.08 ± 0.01 ± 3.0 

1.07 ± 0.01 ± 2.1 

1.08 ± 0.01 ± 2.4 


0.47 

1.05 ± 0.01 ± 4.3 

1.04 ± 0.01 ± 3.6 

1.03 ± 0.01 ± 2.5 

1.03 ± 0.01 ± 0.9 


0.53 

1.01 ± 0.01 ± 3.8 

1.01 ± 0.01 ± 2.3 

1.00 ± 0.01 ± 2.2 

1.03 ± 0.02 ± 1.1 


0.59 

1.00 ± 0.01 ± 5.6 

1.00 ± 0.01 ± 3.7 

0.98 ± 0.01 ± 2.6 

0.99 ± 0.02 ± 1.4 


0.66 

0.96 ± 0.01 ± 5.6 

0.96 ± 0.01 ± 4.7 

0.95 ± 0.01 ± 2.9 

0.95 ± 0.02 ± 3.2 


0.79 

0.94 ± 0.01 ± 5.4 

0.93 ± 0.01 ± 5.5 

0.92 ± 0.01 ± 4.7 

0.89 ± 0.02 ± 1.5 


0.96 

0.92 ± 0.03 ± 10.7 

0.94 ± 0.03 ± 5.0 

0.89 ± 0.02 ± 4.5 

0.85 ± 0.03 ± 3.7 


1.16 

0.82 ± 0.05 ± 7.8 

0.81 ± 0.05 ± 11.9 

0.89 ± 0.04 ± 2.9 

0.81 ± 0.08 ± 10.5 

7T~7r~ 

0.3 

1.07 ± 0.01 ± 3.3 

1.07 ± 0.01 ± 3.2 

1.04 ± 0.01 ± 2.1 

1.02 ± 0.01 ± 1.5 


0.36 

1.11 ± 0.01 ± 2.5 

1.10 ± 0.01 ± 3.2 

1.05 ± 0.01 ± 1.9 

1.07 ± 0.02 ± 2.4 


0.41 

1.09 ± 0.01 ± 3.9 

1.07 ± 0.01 ± 3.1 

1.07 ± 0.01 ± 1.2 

1.08 ± 0.02 ± 1.6 


0.47 

1.02 ± 0.01 ± 4.3 

1.04 ± 0.01 ± 3.2 

1.03 ± 0.01 ± 2.7 

1.04 ± 0.02 ± 2.1 


0.53 

1.01 ± 0.01 ± 4.3 

1.03 ± 0.01 ± 3.4 

0.99 ± 0.01 ± 1.3 

0.99 ± 0.02 ± 1.6 


0.59 

0.99 ± 0.01 ± 5.4 

0.99 ± 0.01 ± 3.9 

0.98 ± 0.01 ± 3.0 

0.99 ± 0.02 ± 4.4 


0.66 

0.96 ± 0.01 ± 5.5 

0.96 ± 0.01 ± 4.9 

0.95 ± 0.01 ± 3.2 

0.96 ± 0.02 ± 1.9 


0.79 

0.91 ± 0.01 ± 6.7 

0.92 ± 0.01 ± 4.3 

0.91 ± 0.01 ± 2.9 

0.91 ± 0.02 ± 2.0 


0.96 

0.86 ± 0.03 ± 5.0 

0.93 ± 0.03 ± 4.3 

0.87 ± 0.02 ± 3.8 

0.85 ± 0.04 ± 3.3 


1.15 

0.73 ± 0.05 ± 8.7 

0.86 ± 0.06 ± 7.6 

0.8 ± 0.05 ± 4.8 

0.79 ± 0.09 ± 6.7 

K+K+ + K-R- 

0.76 

1.07 ± 0.06 ± 5.3 

1.06 ± 0.06 ± 4.5 

1.03 ± 0.05 ± 5.3 

0.98 ± 0.08 ± 7.1 


0.91 

1.05 ± 0.06 ± 8.9 

1.06 ± 0.06 ± 2.9 

1.01 ± 0.05 ± 3.6 

1.05 ± 0.10 ± 6.1 


1.1 

1.08 ± 0.08 ± 5.2 

0.90 ± 0.07 ± 3.8 

0.90 ± 0.07 ± 3.1 

0.98 ± 0.12 ± 2.3 
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TABLE VIII. Azimuthal angle dependence of HBT radii of charged pions, shown as value ± statistical uncertainty [absolute 
value] ± systematic uncertainty [%] for the 0%-20% and 20%-60% centrality bins. 


Centrality 

kr 

[GeV/c] 

0 - 'I'2 

[rad] 

[fm^l 

[fm^] 

Ri^ 

[fm^l 

Ros^ 

[fm^] 

0%-20% 

0.2-0.3 

0 

23.68 ± 0.27 ± 0.56 

25.25 ± 0.34 ± 1.54 

34.72 ± 0.48 ± 1.97 

0.21 ± 0.2 ± 0.11 



7r/4 

22.88 ± 0.27 ± 0.56 

27.33 ± 0.37 ± 1.79 

35.05 ± 0.49 ± 1.63 

2.53 ± 0.21 ± 0.25 



■k/2 

22.46 ± 0.27 ± 0.56 

30.65 ± 0.42 ± 1.8 

35.34 ± 0.51 ± 0.98 

0.26 ± 0.23 ± 0.27 



37r/4 

23.4 ± 0.28 ± 0.51 

28.06 ± 0.39 ± 1.31 

35.46 ± 0.51 ± 1.41 

-2.01 ± 0.22 ± 0.24 


0.3-0.4 

0 

20.5 ± 0.2 ± 0.67 

21.87 ± 0.24 ± 1.38 

24.44 ± 0.29 ± 1.64 

0.22 ± 0.14 ± 0.13 



7r/4 

19.82 ± 0.21 ± 0.4 

24.09 ± 0.27 ± 1.2 

24.03 ± 0.3 ± 1.52 

2.38 ± 0.15 ± 0.26 



7r/2 

18.38 ± 0.2 ± 0.67 

26.93 ± 0.31 ± 2.02 

24.28 ± 0.31 ± 1.38 

0.59 ± 0.16 ± 0.34 



37r/4 

19.16 ± 0.2 ± 0.35 

23.96 ± 0.27 ± 1.78 

24.35 ± 0.3 ± 1.3 

-2.01 ± 0.15 ± 0.31 


0.4-0.5 

0 

17.36 ± 0.18 ± 0.5 

16.61 ± 0.2 ± 1.15 

16.81 ± 0.21 ± 1.31 

0.06 ± 0.12 ± 0.2 



7r/4 

17.01 ± 0.19 ± 0.42 

18.58 ± 0.23 ± 1.41 

17.42 ± 0.23 ± 1.59 

2.11 ± 0.13 ± 0.18 



■KI2 

15.76 ± 0.19 ± 0.32 

20.46 ± 0.27 ± 1.55 

17.73 ± 0.25 ± 1.16 

0.31 ± 0.14 ± 0.19 



37r/4 

16.55 ± 0.19 ± 0.23 

18.51 ± 0.23 ± 1.25 

17.26 ± 0.23 ± 1.29 

-1.79 ± 0.13 ± 0.17 


0.5-0.6 

0 

15.78 ± 0.19 ± 0.45 

13.22 ± 0.19 ± 1.05 

13.85 ± 0.21 ± 0.68 

0.1 ± 0.12 ± 0.25 



7r/4 

15.32 ± 0.2 ± 0.42 

15.34 ± 0.23 ± 1.42 

14.25 ± 0.23 ± 0.88 

1.82 ± 0.13 ± 0.28 



7r/2 

13.75 ± 0.2 ± 0.39 

17.22 ± 0.28 ± 1.31 

14.18 ± 0.25 ± 0.74 

0.16 ± 0.14 ± 0.25 



37r/4 

14.84 ± 0.19 ± 0.43 

14.32 ± 0.22 ± 1.4 

13.85 ± 0.22 ± 0.72 

-1.36 ± 0.12 ± 0.24 


0.6-0.8 

0 

13.87 ± 0.18 ± 0.4 

10.11 ± 0.16 ± 1.04 

9.81 ± 0.16 ± 0.49 

0.04 ± 0.1 ± 0.1 



7r/4 

12.33 ± 0.17 ± 0.35 

11.16 ± 0.18 ± 0.94 

9.88 ± 0.17 ± 1.06 

1.44 ± 0.11 ± 0.26 



7r/2 

11.28 ± 0.17 ± 0.2 

12.01 ± 0.21 ± 1.26 

10.51 ± 0.19 ± 1.24 

0.18 ± 0.11 ± 0.2 



37r/4 

12.35 ± 0.17 ± 0.28 

10.59 ± 0.18 ± 0.78 

9.99 ± 0.17 ± 0.61 

-1.23 ± 0.1 ± 0.17 


0.8-1.5 

0 

11.04 ± 0.24 ± 0.36 

6.85 ± 0.19 ± 0.88 

6.37 ± 0.17 ± 0.74 

-0.07 ± 0.12 ± 0.19 



n j A: 

9.83 ± 0.23 ± 0.31 

7.57 ± 0.23 ± 0.68 

6.05 ± 0.18 ± 0.73 

1.19 ± 0.13 ± 0.15 



7r/2 

7.92 ± 0.21 ± 0.52 

8.61 ± 0.28 ± 0.74 

5.64 ± 0.19 ± 0.44 

0.06 ± 0.13 ± 0.23 



Z-K14, 

9.74 ± 0.25 ± 0.37 

7.04 ± 0.21 ± 0.78 

6.2 ± 0.18 ± 0.48 

-1.28 ± 0.13 ± 0.18 

20%-60% 

0.2-0.3 

0 

16.24 ± 0.19 ± 0.29 

14.53 ± 0.21 ± 0.84 

22.44 ± 0.32 ± 0.97 

0.27 ± 0.13 ± 0.25 



7r/4 

15.18 ± 0.19 ± 0.28 

16.97 ± 0.24 ± 0.85 

22.64 ± 0.33 ± 0.7 

2.3 ± 0.14 ± 0.3 



7r/2 

12.81 ± 0.17 ± 0.24 

18.97 ± 0.28 ± 0.88 

21.85 ± 0.34 ± 0.76 

0.3 ± 0.14 ± 0.34 



3'7r/4 

14.79 ± 0.18 ± 0.33 

16.34 ± 0.24 ± 0.95 

22.47 ± 0.33 ± 0.52 

-1.8 ± 0.13 ± 0.35 


0.3-0.4 

0 

13.98 ± 0.15 ± 0.27 

12.92 ± 0.14 ± 0.62 

15.9 ± 0.2 ± 0.59 

0.09 ± 0.09 ± 0.15 



7r/4 

13.15 ± 0.15 ± 0.22 

15.38 ± 0.18 ± 0.79 

15.77 ± 0.21 ± 0.46 

2.82 ± 0.1 ± 0.12 



7^12 

11.04 ± 0.13 ± 0.26 

17.49 ± 0.22 ± 0.74 

15.11 ± 0.21 ± 0.52 

0.31 ± 0.11 ± 0.15 



37r/4 

12.56 ± 0.14 ± 0.18 

15.03 ± 0.17 ± 0.74 

15.19 ± 0.2 ± 0.63 

-2.4 ± 0.1 ± 0.16 


0.4-0.5 

0 

12.58 ± 0.14 ± 0.17 

10.22 ± 0.12 ± 0.58 

11.31 ± 0.15 ± 0.45 

0.17 ± 0.08 ± 0.18 



7r/4 

11.55 ± 0.14 ± 0.22 

12.12 ± 0.15 ± 0.64 

11.31 ± 0.16 ± 0.51 

2.49 ± 0.09 ± 0.14 



7r/2 

9.32 ± 0.13 ± 0.28 

14.99 ± 0.21 ± 0.82 

11.39 ± 0.18 ± 0.63 

0.17 ± 0.1 ± 0.13 



37r/4 

11.05 ± 0.13 ± 0.26 

12.25 ± 0.16 ± 0.58 

11.72 ± 0.17 ± 0.59 

-2.04 ± 0.09 ± 0.19 


0.5-0.6 

0 

11.28 ± 0.14 ± 0.17 

8.22 ± 0.12 ± 0.47 

9.08 ± 0.14 ± 0.29 

0.09 ± 0.09 ± 0.1 



7r/4 

10.15 ± 0.14 ± 0.29 

9.83 ± 0.15 ± 0.48 

9.38 ± 0.16 ± 0.39 

2.21 ± 0.09 ± 0.18 



7r/2 

8.12 ± 0.14 ± 0.27 

12.04 ± 0.21 ± 0.56 

8.96 ± 0.18 ± 0.49 

-0.03 ± 0.1 ± 0.16 



37r/4 

10.09 ± 0.14 ± 0.24 

10.2 ± 0.16 ± 0.64 

9.07 ± 0.16 ± 0.41 

-2.05 ± 0.09 ± 0.13 


0.6-0.8 

0 

9.57 ± 0.12 ± 0.21 

6.08 ± 0.09 ± 0.4 

6.53 ± 0.1 ± 0.31 

0 ± 0.06 ± 0.06 



7r/4 

8.55 ± 0.13 ± 0.12 

7.68 ± 0.13 ± 0.41 

6.46 ± 0.12 ± 0.25 

1.82 ± 0.07 ± 0.12 



7r/2 

6.71 ± 0.12 ± 0.31 

9.47 ± 0.18 ± 0.47 

6.61 ± 0.14 ± 0.27 

0.16 ± 0.08 ± 0.14 



37r/4 

8.62 ± 0.13 ± 0.26 

7.56 ± 0.13 ± 0.45 

6.57 ± 0.12 ± 0.44 

-1.82 ± 0.07 ± 0.16 


0.8-1.5 

0 

7.81 ± 0.16 ± 0.27 

4.2 ± 0.1 ± 0.29 

4 ± 0.1 ± 0.2 

0 ± 0.07 ± 0.1 



7r/4 

7.18 ± 0.18 ± 0.33 

5.06 ± 0.15 ± 0.27 

4.3 ± 0.13 ± 0.24 

1.47 ± 0.09 ± 0.12 



7r/2 

4.74 ± 0.16 ± 0.22 

6.72 ± 0.26 ± 0.68 

3.84 ± 0.15 ± 0.65 

0.34 ± 0.1 ± 0.2 



37r/4 

6.58 ± 0.17 ± 0.36 

5.36 ± 0.16 ± 0.25 

3.94 ± 0.12 ± 0.45 

-1.32 ± 0.09 ± 0.09 
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TABLE IX. Azimuthal angle dependence of HBT radii of charged kaons, shown as value ± statistical uncertainty [absolute 
value] ± systematic uncertainty [%] for the centrality bins shown in Fig. (8] 


Centrality 

(/) - ^<2 

Rs^ 

Ro^ 

Ri^ 

Ros^ 


[rad] 

[fm^] 

[fm^] 

[fm^] 

[fm^] 

0%-20% 

0 

11.86 ± 0.68 ± 0.42 

10.77 ± 0.79 ± 1.14 

10.43 ± 0.72 ± 1.34 

1.22 ± 0.58 ± 0.77 


7r/4 

10.02 ± 0.57 ± 1.04 

11.89 ± 0.82 ± 1.22 

8.74 ± 0.61 ± 1.70 

2.43 ± 0.52 ± 0.65 


7r/2 

7.98 ± 0.47 ± 0.92 

13.29 ± 0.97 ± 1.96 

9.25 ± 0.64 ± 1.47 

-1.31 ± 0.52 ± 0.42 


37r/4 

10.67 ± 0.61 ± 0.66 

12.45 ± 0.90 ± 0.86 

9.44 ± 0.65 ± 1.09 

-2.15 ± 0.61 ± -0.69 

20%-60% 

0 

10 ± 0.58 ± 0.69 

7.49 ± 0.54 ± 0.69 

5.71 ± 0.40 ± 0.6 

0.28 ± 0.44 ± 0.22 


7r/4 

7.52 ± 0.46 ± 0.56 

8.03 ± 0.58 ± 0.69 

6.07 ± 0.43 ± 0.71 

2.05 ± 0.40 ± 0.29 


7r/2 

4.69 ± 0.31 ± 0.68 

8.56 ± 0.59 ± 0.99 

5.61 ± 0.42 ± 0.81 

0.41 ± 0.34 ± 0.25 


37r/4 

8.03 ± 0.49 ± 0.65 

8.51 ± 0.62 ± 0.77 

5.27 ± 0.40 ± 0.83 

-1.52 ± 0.43 ± -0.24 


TABLE X. Oscillation amplitudes relative to the event plane for charged pions, shown 
[absolute value] ± systematic uncertainty [%] for the 0%-20% and 20%-60% centrality bins 

as 

value 

± statistical uncertainty 

Centrality 

rriT 

[GeV/c] 



[fm^] 





Ro^ 

[fm^] 





Ri^ 

[fm^l 





Ros^ 

[fm^] 



0%-20% 

0.30 

0.026 

± 

0.007 

± 

0.015 

0.095 

± 

0.008 

± 

0.024 

0.114 

± 

0.010 

± 

0.029 

0.099 

± 

0.007 

± 

0.009 


0.38 

0.054 

± 

0.006 

± 

0.032 

0.101 

± 

0.007 

± 

0.018 

0.125 

± 

0.008 

± 

0.025 

0.113 

± 

0.005 

± 

0.013 


0.47 

0.046 

± 

0.007 

± 

0.018 

0.104 

± 

0.008 

± 

0.025 

0.116 

± 

0.008 

± 

0.030 

0.117 

± 

0.006 

± 

0.010 


0.57 

0.065 

± 

0.008 

± 

0.027 

0.125 

± 

0.009 

± 

0.020 

0.126 

± 

0.009 

± 

0.024 

0.107 

± 

0.006 

± 

0.016 


0.70 

0.106 

± 

0.008 

± 

0.015 

0.082 

± 

0.010 

± 

0.032 

0.072 

± 

0.009 

± 

0.028 

0.107 

± 

0.006 

± 

0.016 


0.93 

0.160 

± 

0.014 

± 

0.040 

0.098 

± 

0.019 

± 

0.040 

0.077 

± 

0.015 

± 

0.022 

0.129 

± 

0.010 

± 

0.013 

20%-60% 

0.30 

0.114 

± 

0.007 

± 

0.009 

0.131 

± 

0.009 

± 

0.024 

0.149 

± 

0.010 

± 

0.027 

0.140 

± 

0.007 

± 

0.022 


0.38 

0.113 

± 

0.007 

± 

0.019 

0.150 

± 

0.007 

± 

0.015 

0.181 

± 

0.008 

± 

0.020 

0.207 

± 

0.006 

± 

0.010 


0.47 

0.144 

± 

0.007 

± 

0.017 

0.180 

± 

0.008 

± 

0.015 

0.201 

± 

0.009 

± 

0.019 

0.205 

± 

0.006 

± 

0.012 


0.57 

0.155 

± 

0.008 

± 

0.022 

0.185 

± 

0.010 

± 

0.011 

0.189 

± 

0.010 

± 

0.011 

0.216 

± 

0.007 

± 

0.014 


0.70 

0.165 

± 

0.009 

± 

0.029 

0.212 

± 

0.010 

± 

0.021 

0.196 

± 

0.010 

± 

0.018 

0.220 

± 

0.007 

± 

0.017 


0.93 

0.226 

± 

0.015 

± 

0.026 

0.211 

± 

0.019 

± 

0.027 

0.172 

± 

0.016 

± 

0.021 

0.216 

± 

0.010 

± 

0.012 


TABLE XL Oscillation amplitudes relative to the event plane for charged kaons, ± systematic uncertainty [%] for the 0%-20% 
and 20%-60% centrality bins shown in Fig. [T^j 


Centrality 

rriT 

R.^ 



p 2 
-^os 


[GeV/c] 

[fm^] 

[fm^] 

[fm^l 

[fm^l 

0%-20% 

0.91 

0.193 ± 0.034 ± 0.078 

0.106 ± 0.044 ± 0.033 

0.128 ± 0.052 ± 0.044 

0.227 ± 0.040 ± 0.067 

20%-60% 

0.91 

0.360 ± 0.035 ± 0.080 

0.070 ± 0.042 ± 0.035 

0.076 ± 0.045 ± 0.038 

0.238 ± 0.040 ± 0.030 
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